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Abstract. We establish a link between the dual parametrization of GPDs and a popular parametrization 
based on the double distribution Ansatz, which is in prevalent use in phenomenological applications. 
We compute several first forward-like functions that express the double distribution Ansatz for GPDs in 
the framework of the dual parametrization and show that these forward-like functions make the dominant 
contribution into the GPD quintessence function. We also argue that the forward-like functions Q2v (x) with 
v > 1 contribute to the leading singular small- a; behavior of the imaginary part of DVCS amplitude. This 
makes the small-XBj behavior of ImA^^^^ independent of the asymptotic behavior of PDFs. Assuming 
analyticity of Mellin moments of GPDs in the Mellin space we are able to fix the value of the D-form 
factor in terms of the GPD quintessence function N{x,t) and the forward-like function Qo{x,t). 

PACS. 13.60.Fz Elastic and Compton scattering - 11.55.Hx Sum rules 



1 Introduction 

Generalized parton distributions (GPDs) [1] have been in- 
vestigated in great details during the past decade. These 
distributions, proved to be extremely efhcient for the de- 
scription of the quark and gluon structure of hadrons. 
GPDs provide a natural generalization of parton distri- 
bution functions, familiar from inclusive reactions, and 
elastic form-factors. The current understanding of both 
theoretical and experimental aspects of GPDs is reviewed 
in refs. [2l[3l|4l[5] . Extraction of GPDs from the experi- 
mental data is highly demanded, since these functions 
contain valuable new information on hadron structure. 
The theoretical opportunity to access GPDs experimen- 
tally is provided by the coUinear factorization theorems for 
hard exclusive reactions [6l[7]. The dedicated experiments 
[8] performed during the last years provide an increas- 
ing amount of precise experimental data. Unfortunately, 
since GPDs are complicated functions of the longitudinal 
momentum fraction of partons (x), skewness parameter 
(^), the momentum transfer squared (t) as well as of fac- 
torization scale, the direct extraction of GPDs from the 
observables turns out to be an extremely difficult task. 
Moreover, GPDs always enter the observable quantities 
being integrated over x with weighting functions given by 
the propagators of partons between the incoming virtual 



Send offprint requests to: 
" Maxim.Polyakov@tp2.ruhr-uni-bochum.de 
semenovk@tp2.ruhr-uni-bochum.de 



photon and the outgoing real photon or meson. There- 
fore, in order to extract GPDs from the data, one usually 
rely on different phenomenologically motivated parame- 
terizations and simultaneous fitting procedures for several 
observables. 

An important advantage of the dual parametrization 
of GPDs [9J is that it suggests the possible form of de- 
convolution procedure [IOlfTT[|T2] which allows to specify 
the maximum amount of information on GPDs that can 
be extracted from the experimental data in a single GPD 
quintessence function. GPD quintessence function can be 
unambiguously recovered from the leading order ampli- 
tude with the help of the Abel transform tomography 
method. Another gain from the dual parametrization is 
that it allows one to distinguish the contribution to the 
observables brought by parton densities from the genuine 
non- forward contributions (see discussion in [13j). 

In this paper, using reparametrization procedure de- 
scribed in [T^, we establish the link between the dual 
parametrization of GPDs and the famous Radyushkin dou- 
ble distribution Ansatz [l5lfT6lfT7[|18[fT9] employed in nu- 
merous phenomenological applications. The explicit ex- 
pressions for the first forward-like functions Q2 and Q4 al- 
low us to quantify the relative importance of non-forward 
effects encoded in these functions. We also argue that, in 
the framework of the dual parametrization, the forward- 
like functions Q2u{x) with z/ > 1 may contribute to the 
leading singular small-^ asymptotic behavior of the imag- 
inary part of the leading order DVCS amplitude. In this 
case the forward-like functions Q2u{x) with v > 1 may 
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have the small-x singularities, which lead to divergencies 
of the generaUzed form- factors i?2y-i o- We discuss the 
regularization procedure that allows to assign a precise 
meaning to the potentially divergent integrals for gener- 
alized form factors B2U-1 o- Assuming the analyticity of 
Mellin moments of GPD in Mellin space (see [^[2T[|22] ) 
one can unambiguously fix the contribution of the _D-form 
factor into the real part of the DVCS amplitude in terms 
of the GPD quintessence function and the forward-like 
function Qo- 

Finally, in the Appendix [X] in order to avoid confu- 
sions in the Hterature we review the form of the inte- 
gral transformation [9lfT0] relating GPD and the set of 
forward-like functions Q2U in the framework of the dual 
parametrization. We generalize this result for the case 
when the forward- like functions Q2v{x) with v > 1 are 
allowed to have small- a; singularities, which lead to diver- 
gencies of the generalized form- factors o- 



2 Basic facts on the dual parametrization of 
GPDs 

The dual parametrization of GPDs [S] is based on the 
representation of GPDs as infinite sums of t-channel reso- 
nance exchanges j23j . Originally, the dual parametrization 
was formulated for the case of spinless hadrons. First we 
briefly discuss the difference, which arises when dealing 
with spin one half particles (a detailed discussion on this 
issue is presented in [Mj)- Note that the specifics related to 
spin-i was not taken into account in the early phenomeno- 
logical applications of the dual parametrization [25 , as it 
has been stressed in [26] . 

According to [51l27j. the following (so-called electric 
and magnetic) combinations of nuclcon GPDs H'^{x,^,t) 
and E'^(x,£,,t) are suitable for partial wave expansion in 
the t-channel: 



t 



(x, e, t) = F?(x, t) + -r^E'^ix, C, t) , 



4:171% 



H^^'\x,^,t) ^ H'^ix,i,t) + E-^ix^tt) ■ 



(1) 



Here we employ the nucleon generalized parton distribu- 
tion of the particular flavor q: H''{x,^,t), E'^{x,^,t) with 
the properties listed below. 

— Generahzed parton distributions H''{x, ^, t), E'^{x, t) 
defined for x S [—1,1] are reduced to the following t- 
dependent quark densities in the limit ^ 0: 



iJ'?(x,e = 0,i) 



ii;«(x,C = 0,t) 



q(x, t) for X > 
—q{—x, t) for a; < 



e'?(x,i) for x>0 
-e«(-x,t) for a; < 



(2) 



(3) 



— The following normalization conditions are fulfilled: 

/•I 4 
/ dxxH^^'>{x,^,t^Q)^Ml + -dU^; 
J-i 5 

J dxxH^^\x,^,t = 0) = 2J\ (4) 

where Ml stands for the momentum fraction carried 
by quarks and antiquarks of the particular flavor q in 
the nucleon; J* denotes angular momentum carried 
by quarks of flavor q; d\ is the first coefficient of the 
Gegenbauer expansion of the _D-tcrm of fiavor q. 

The dual parametrization can be introduced for the 
singlet and nonsinglet combinations of Below we 

consider only the case of the singlet (C = +1) combina- 
tions of H^^'^\ which enter the description of DVCS. 
These combinations are introduced according to 

h(''\x,u)-h^''\~x,U); 

--H("'\x,^,t)-H('^'\-x,tt). (5) 



Hr>{x,^,t) 



Hr'{x,tt) 



Note that the singlet electric and magnetic combinations 
of nucleon GPDs iJ (^.^^) (x, t) are defined for x e [0, 1] . 
For ^ = they are reduced to the following combinations 
of ^-dependent quark densities q+{x,t) = q{x,t) + q{x,t), 
e\{x,t) = e9(x,t) -|-e9(x,t): 



Hf\x,^^Q,t)^q+{x,t) + 



H':^'\x,i^O,t)^q+{x,t)+e%{x,t). 



(6) 
(7) 



The singlet electric and magnetic combinations of nucleon 
GPDs n'f''^^ {x, ^, t) are normalized according to 



dxxHf\x,S,,t^ 0) 



Ml 



-die- 

5 



dxa;i7f^^(2;,C,t = 0) = 2J«. 



(8) 



All formulae for the electric combination of nucleon 
GPDs are the same as in the case of the scalar hadrons 
addressed in [9j. Thus, the partial wave decomposition in 
the t-channel for the singlet electric combination of GPDs 
H^\x,^,t) is written as the following formal seriefl: 



Hi''\x,U)^ 

OQ n+1 

Tl=l 1=0 

odd even 



1-^ 

e 



Pi ; 



(9) 



^ Note that we have introduced an additional factor 2 in the 
partial wave expansions (|10p . This is done to make the re- 
sulting GPDs Hf''^''' satisfy (0), © in the limit ^ 0. In the 
original paper ^9 a rather unusual convention that differs by 
a factor i from (O was employed. In particular this unusual 
convention has led to much confusion in the early phenomeno- 
logical applications of the dual parametrization [281125] (see [29] 
for the detailed discussion of this issue). 
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where Cn (x) stand for the Gegenbauer polynomials; Pi (x) 
are Legendre polynomials; Bni{t) are the generalized form 
factors; x, ^ and t stand for usual GPD variables. As it is 
pointed out in [24j . the i-channel partial wave expansion 

for the magnetic combination goes over |-P/ 



and by the _D-form factor D"^, that is given by 



oo n+1 



X 

e 



ri2 



az 

-1 



oo 



1 



Z 



where 



{l-z^)Y^dlCl'\z) 



(15) 



(16) 



1 

1 

(10) 



In the following we consider only the case of the singlet 
(C = +1) electric nucleon GPD. The variable t plays no 
particular role for our analysis, so for simplicity we set 
it to zero for the rest of this section. In this limit the 
singlet electric combination of nucleon GPDs is reduced 
to the usual C = +1 nucleon GPD Thus, we omit the 
superscript "(£')": 

H^f{x,^,t^Q) = H+{x,0- 

In the forward limit H^{x,S,) is reduced to q^(x), where 
q+{x) = q{x) + q{x) stands for the singlet combination 
of forward quark distributions. However, the generaliza- 
tion of our analysis (in particular, of the main results of 
sects. [3l HJ for fixed t ^ is straightforward. 

Below, we list a number of important facts concerning 
the dual parametrization of GPDs. 

— In order to sum up the formal series ^ for _ff+(x,^) 
(see Appendix]^ one has to introduce the set of forward- 
like functions Q2v{x) (y = 0, 1, ...), whose Mellin mo- 
ments generate the generalized form factors BnV- 



B 



n n+1 — 2i/ 



dxx'^Q2u{x) . 



(11) 



— At the LO the scale dependence of the forward-like 
functions Q2i/{x) is given by the standard DGLAP evo- 
lution equation so that these functions evolve as usual 
parton distributions under QCD evolution. 

— The forward-like function Qo{x) is related to the for- 
ward quark densities according to 

Q,{x)=q+{x)-^ f %q+{y). (12) 

2 Jx y 

— The leading order twist-2 amplitude A{^) = A{£_,t — 
0): 



dx H+{x, ^) 



(13) 

is completely determined by the GPD-quintessence func- 
tion N{x) [lOlfTT]: 



N{x) 



v=0 



(14) 



stands for the _D-term. The calculations performed in 
the chiral quark-soliton model [30] suggest D^{z) « 
D'^{z), thus D^^iz) = ■^D{z). Nf stands here for the 

number of flavors and D{z) — (1— z^) J2n odd dnCn^^{z 
is the flavor singlet D-term. 

The partial wave expansion of the DVCS amplitude 
([T51) reads [S] 



oo n+1 

n = l 1 = 



1 

n iPl I ^ 



(17) 



This formal series can be summed up yielding the fol- 
lowing result for the DVCS amplitude: 



A{0 



J OO 

-E^(^ 



1 



1 

/■^ dx 

+4 / —{N{x)~Qo{x)). 
Jo X 



2x 



+ 



2x 



(18) 



The explicit expressions for the imaginary and real 
parts of the DVCS amplitude ^ read [9lfT0! 



lmA{0 = 2 



dx 



N{x) 



1 



(19) 



ReA{0 = 4£i« -f 2 



dx 



N{x) 



2x_ 



1 + ^ 



dx , 
— N{x) 

X 



1 + ^ 



(20) 



The ZJ-form factor D"^ can be formalljQ expressed through 
the GPD quintessence function N{x) and the forward- 



The problem of treating the possible divergencies in (|2ip is 
discussed in sect. [3] 
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like function Qq according to [10] 



—N{x) 



+ 



dx 



(21) 



{N{x)-Qo{x)) 



The polynomiality condition for Mellin moments of 
generalized parton distribution Ij implies that for N = 
1,3,...: 



[ dxx^H+{x,C) 
Jo 



jAN) fN+1 



(22) 

In the framework of the dual parametrization, the set 
of coefficients /ig^' with i/ < -^^^ is given by 



AN) 



N n+1 

EE^"' (-1)' 

11 = 1 i = 



(n+l) (n + 2) r{N+l) 



N-r 



(23) 



Note that due to the Gamma function in the denom- 
inator only the terms with / > N + 1 — 2iy are non 
zero in the sum over I in (|23p . As the result for the 
given odd N and i' {i' < ^^y^) only the forward- like 
functions Q2p, with /i = 0, 1, ... ,1^ are relevant for the 

computation of the coefficient h^K 
An important role is played by the expansion of GPD 
iJ+(a;,^) in powers of ^ around the point C = with 
fixed a; (x > ^). To the order this expansion is given 

by m 

H+{x,0 = Hf\x)+eHi'\x)+eH^+\x) + ... 
dy 



Qo{x) 



yx 
2 

x^ d 



y 



3/2 



Qo{y) 



Ax dx 



1 

32 



dyQoiy) 



yi 

_3_ 

32 



dyQ2{y)- 



of forward-like functions through GPDs for various phe- 
nomenological parametrizations of GPDs. This enables us 
to rewrite the particular phenomenological parametriza- 
tions of GPDs in the framework of the dual parametriza- 
tion [IJ. Let us assume that the expansion of GPD H+(x, ^) 
around ^ = for a: > ^ calculated in the framework of a 
certain model is known: 



H+{x,o = Mx) + Mx)e + Mx)e + oif) . 

where 



(25) 



1 



-^H+{x,C)i=Q- 



{2v)\ dC 

Below we list the explicit expressions for several first forward- 
like functions Q2 {x) and Q4 (x) derived in |14j from match- 
ing the expansions ([M)) and (P5)) . Clearly, since the GPD 
calculated in the realistic model is supposed to have the 
correct forward limit 

(t)Q{x) EE H+{x, ^ = 0) = q+{x) , 

the usual result (fT2ll for Qo{x) is recovered. For Q2{x) 
such expression reads 



Q2{x) 
r-1 

+ ' 



2{l~x^) 



dy 

+4.Mx) 



x^ 
-Ibx 



q+ix) 



(l-x^) , 



4y4 



3 5 X . 



dyMy) 



15x 
4y2 



J3_ 
2^ 



3_ 

4x 



(26) 



The explicit expression for Q4{x) reads 

Qiix) ■ 



49 
8x^ 



11 

4 a; 



7 

4x3 
1 

dy q+{y) 



17 

8x2 

(1-. 

/315x 
\16y^ 



(1 - x^) g+(x) 

(l-x2) 



+ [^-20]Mx) 



35 



2x2 

135 X 21 
8 2/4 ~ 4^ 



qUx) 

27 X 



16 y^ 



(24) 



21 
4^ 



dy 02 (y) 

15 

~ 8x 



X 

7 



dy4>i{y) z~ 



315x 

16/ ~ 16x3 

4>4ix) 

315x 



02 (x) 

35 
4y3 



135 X 



3 

45 

X 



16 y2 



35 



15y ^ 35 y2 
4x2 16x3 



(27) 

The result up to the order was derived in [13]. A re- 
markable property of such expansion is that up to the In principle, the derivation of the analogous expressions 
particular order it involves only a finite number of for forward-like functions Q21, with v > 3 is straightfor- 
forward-Hke functions Q2ti{x) with jj, < v. This property ward. Unfortunately, the corresponding expressions turn 
allows to invert the expansion (|24p and to express the set out to be too bulky. 
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3 D-term and the divergencies of generalized 
form factors 

In the framework of the dual parametrization of GPDs 
the small-^ behavior of the imaginary part of DVCS am- 
plitude (fT9|) inherits the singular behavior of the GPD 
quintessence function N{x) for small x. Assuming that 
N{x) ^ for a; ~ 0, one can establish the following 
asymptotic behavior of the imaginary part of the leading 
order DVCS amplitude ^ [9]: 



W(0 



2"+ir(i)r(a + i) 

1) 



(28) 



It is usually supposed, that the index a of the leading 
singular power-like asymptotic behavior of the imaginary 
part of DVCS amplitude for small ^ is determined by 
the singular behavior of the corresponding forward sin- 
glet quark distributions measured in DIS experiments: 
q+{x) ~ for X (0 < a < 2, a ^ 1). In the 
phenomenological applications of the dual parametriza- 
tion pSlf^ it was tacitly assumed that the the leading 
small-^ singular behavior of the imaginary part of DVCS 
amplitude is entirely determined by the small- a; behavior 
of the forward- like function Qo{x): 



lix) 

x" 



1 



a -I- 1 x° 



(29) 



The functions Q2u{x) with high index v were supposed 
to be appropriately suppressed for small x thus making 
no influence on the small-^ singular behavior of Imyl(^). 
However, the corresponding assumption is too restrictive 
and by no means motivated from the physical point of 
view. The forward- like functions Q2u{x) with v > 1, in 
principle, may contribute to the coefflcient at the leading 
singular power of ^ in the asymptotic expansion of ImA(^) 
for small ^. For this one has to assume the following lead- 
ing small- a: behavior of Q2v{x) with v >1: 



2i^{x) 



1 



(30) 



where 0<a<2(Q;7^1).A problem that immediately 
arises in this case is that the asymptotic behavior (|30p 
of Q2u{x) with V > I leads to the divergencies of the 
generalized form factors i?2i/-i 



f-x^^Q2.[x). 
Jo ^ 



2u-l — 



(31) 



According to pTjl . the generalized form factors PT|) are 
the lowest order Mellin moments of the forward-like func- 
tions Q2v with V > Q relevant for the calculation of GPDs 
in the framework of the dual parametrization. The cru- 
cial point is that as it can be easily seen from (fT8| , in the 
calculation of the leading order DVCS amplitude the po- 
tentially divergent generalized form factors (|3ip contribute 
only into the _D-form factor (|2ip . 



^ Below for simplicity we set the variable t to zero. The gen- 
eralization of our analysis for t 7^ is straightforward. 



Thus, we have to specify a suitable regularization of a 
possible non-integrable singularity at a; = in (|3ip . We 
restrict our following analysis to the case in which the 
contributions x^'^Q2u{x) of the forward-like functions into 
GPD quintessence function A''(a;) as well as A''(a;) itself 
belong to the class of functions with power like behavior 
for a; '--^ 0, which can be presented as the finite sums of 
singular terms [31] : 



1 



2. (a;), N{x) elF: F{x) = ^ —fr{x) 



r=l 



(32) 



Here frix) with r = 1, ...R stand for arbitrary functions of 
X that are infinitely differentiable in the vicinity of a; = 0. 
It is also supposed that frix) have zeroes of a sufficiently 
high order for a; = 1. We assume that < < 2 (a^ ^ 1) 
so that a non integrable singularity at x = manifests it- 
self only in the computation of the lowest order Mellin 
moments ((3T|) of Q2iy{x) [v > 1) relevant for the calcula- 
tion of GPD in the dual parametrization. All higher order 
Mellin moments Bn n+i-2L' (fTTjl with n = 2i'+l, 2z/-|-3, ... 
turn to be finite. The class of functions defined in Eq. (p2)) 
seems to be convenient for model building. In particular, 
the functions usually employed in fitting of PDFs belong 
to this class. Moreover, the imaginary part of the DVCS 
amplitude ImA(^) computed using the GPD quintessence 
function N{x) from the class l\'3'2\i also belongs to this class. 

In order to work out a method of handling the diver- 
gencies of the generalized form factors ([3T|) one has to 
consider the general analytic properties of the DVCS am- 
plitude. According to the analysis presented in '321 l331l?fl 
IM] the leading order DVCS amplitude is considered to be 
a holomorphic function of the variable uj ~ ^. The real 
part of the amplitude can be expressed through its imag- 
inary part with the help of a single variable dispersion 
relation in ui for fixed value of t. In the singlet case this 
dispersion relation require one subtraction. The subtrac- 
tion constant is given by the value of the amplitude at 
the non-physical point w = (^ = 00). It is known to 
be fixed by the D term and equals where stands 
for the D-form factor (fT5|) . After one switches back to the 
variable f , the once subtracted dispersion relation for the 
DVCS amplitude ^ reads 



A(e) = 4L>« 



$ - e - «e e + e - *e 



ImA(C' - ie) 



(33) 



In general, the subtraction constant in a dispersion rela- 
tion presents an independent quantity, which cannot be 
fixed just with help of the information on the disconti- 
nuities of the amplitude. Thus, in order to determine the 
value of the subtraction constant one has to attain certain 
additional information on the amplitude under considera- 
tion. 

The problem of fixing the subtraction constants in the 
dispersion relation for the amplitude solely in terms of 
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the absorptive part of the amplitude is addressed in the 
theory of the analytic 5-niatrix. The key for this issue is 
provided by assuming the so-called analyticity of the sec- 
ond kind that is the analyticity of the amplitude in the 
complex angular momentum plane [35j . The consideration 
of the analytical properties of the scattering amplitude in 
the complex angular momentum plane is usually associ- 
ated with the partial wave expansion of the amplitude 
[36] , However, as pointed out in [37] (see also discussion 
in [35]), this type of analysis can equally be applied for the 
expansions of the amplitude other than the partial wave 
expansion in particular for the expansion of the amplitude 
into the ordinary power series. 

This very logics was employed in [20l[2Tll22] where it 
was suggested to fix the subtraction constant in the dis- 
persion relation (j33p in terms of the imaginary part of the 
DVCS amplitude assuming analyticity properties in j of 
the specific combinations of coefficients 'bS^^^^^ at powers 
of ^ of the Mellin moments of GPD (see eq. ([^^ for the 
definition) . 

For this following the line of analysis of [20 l l2i p 22j let 
us consider the family of GPD sum rules for the coef- 
ficients at powers of ^ of the Mellin moments of GPD 

together 
32,33,27, 



stant in the dispersion relation (|33p [21] : 



'"21/ 

with 



Employing the dispersion relation 
131 one can establish the GPD sum rule 



dx 



1 



1 



[H+{x,0-H+{x,x)]=^D'' . 



(34) 



Here the principle value prescription is dropped since the 
singularity at a; = ^ turns to be integrable. Expanding the 
convolution kernel in (|34]l in powers of j and employing 
the polynomiality property (|22)) of the Mellin moments of 
GPD one derives a family of sum rules for the coefficients 
at powers of ^ of Mellin moments of GPD: 



oo „i 

/ dxx^[H+{x, 



a;)-H+(a;,0)] , 



(35) 



with j = 1, 3, .... Note that under the usual assumptions 
on the small- a: asymptotic behavior: 

i?+(a:,a:) - i7+(a;,0) - 4^ with a<2 ia^X) 

(36) 

the integral in (|35]) actually converges for odd positive j. 
The substantial step is to define the function 

#(j) = / dxx^ [H+{x,x) - H+{x,0)] 



(38) 



In order to take advantage of the analyticity of ^>(j) in the 
complex j plane we have to specify our assumptions for the 
analytic properties of the function H+{x, x) — H+{x, 0) = 
i [ImA(2;) — IuyAdis{x)\ in the x plane. We suppose the 
Regge-like asymptotic behavior ([M]l oilT[iA{x)~lv[\ADi si^) 
for small- a; which can be established experimentally. In 
our present analysis in order to obtain simple analytical 
properties of <?(j) in the complex j plane we assume that 
lxa.A{x) — ImADisix) belongs to the class of functions 
defined in Eq. (j32[j i. Under this assumption the func- 
tion <?(j) turns to be a meromorphic function of j for 
Rcj > —1 — e (e > 0). <?(j) has just a finite number of 
simple poles in the complex plane j for Rej > —1. In this 
case the analytic continuation of <?(j) to j = —1 exists 
and is obviously unique. It is provided by the following 
explicit expression for the regularization of the divergent 
integrals in ([37| : 



Reg tdxi^.tdxi^ 

' dx^[fix)^f{0)~xf'm 







X 



-/(0)-^/'(0)^ 



(39) 



Clearly this is what is known as the analytic regularization 
of integrals [31] . 

Thus, under the specified above analyticity assump- 
tions (which, however, can be violated, as we discuss be- 
low), the subtraction constant in the dispersion relation 
(|33p can be fixed in terms of the imaginary part of the 
DVCS amplitude according to [12] 



1 

201^ dx-\H+ix,x) - H+(x,0)] . 
J(o) X 



(40) 



Now we discuss how these general considerations ap- 
ply for the case of the dual parametrization of GPDs. We 
argue that the implementation of analyticity of the combi- 
nations (|34)) of coefficients of the Mellin moments of GPDs 
suggests the use of analytic regularization for the poten- 
tially divergent generalized form factors i32i/^i o (|?T]) . 

Let us consider the sum rule (|40)) for the D- form factor 
in the framework of the dual parametrization of GPDs. 
The explicit expressions for H{x,x) and H{x,0) through 
the forward like functions in the framework of the dual 
parametrization read 



1 /-i" . H{x,0) = Qo{x) + ^ [' ^Qo{y) 

- dxx^ [lmA{x) ^ImADisix)] . (37) 2 y2 

TT Jo 



(41) 



Once the function <?(j) is analytically continued to j = —1 
it provides the desired relation for the subtraction con- 



* The generalization for ImA(x) — lmADis{x), which be- 
long to the more intricate classes like e.g. F : F(x) = 
-^{logxY^ fr{x) is straightforward. 
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and 



H{x,x =- / ^ 



N{y) 



2y _y2 



(42) 



1 



We assume that both functions N{y) and Qo{y) belong 
to the class (|32|) . The leading singular behavior for y '--^ 
is Qo{y), N{y) ^ ^ with a < 2. It is straightforward to 
check that interchanging the order of integration 



2 dx 



7i'J(o) y JO 



V 



i+y' dx 



22/ , 

X 

1 



2y 



(43) 



y 



1 



is a rigorously defined operation for the class of functions 
under consideration. The same is true for 



(0) 



dx yjx 

~x~ 



' dy. 



My) 



X y- 



1 dy^ 

(0) yi 



My) 



dx 



(44) 



In this way one arrives to the analytically regularized ver- 
sion of the expression (PTjl for the D form factoiQ: 



y '"'\ViT^ 

' ^ [N{y) - Qo{y)] . 

(0) y 



(45) 



The integral in the first term of (|45p converges under our 
assumptions on the small-y behavior of N(y) so it does 
not require analytic regularization. The second term is the 
pure contribution of the problematic form factors B2U-1 
since 



1 dy 

(0) y 



N{y) 



Qo{y)]^ t dyf2-y'-Q2Ay). 



So we conclude that the analytic regularization of the in- 
tegral in the second term of (|45l) suggests the use of the 
analytic regularization for the individual terms of the sum 



= / —x'^''Q2uix) ■ 
J{0) X 



(46) 



It is also extremely instructive to check that under our 
assumption on the analyticity of the combinations (|37|) of 
coefficients of the Mellin moments GPDs in the framework 
of the dual parametrization satisfy the so-called "dual- 
ity property" suggested in [22] . According to the "duality 
property" , GPD in the central region < x < ^ can be 



^ In [TT] an alternative regularization prescription was sug- 
gested. The two regularizations obviously differ just by a finite 
constant and thus are equivalent up to a finite _D-term contri- 
bution. 



completely restored from its knowledge in the outer re- 
gion ^ < a; < 1. The coefficients at powers of ^ of N-th 
Mellin moment {N = 1, 3, ...) of GPD H+{x,l) can be 
determined from the small ^ expansion (|24| of H^(x,^) 
in powers of ^ for a; > ^ according to 



j/=a 

t^r, Jo J(Q) 



10 

N-1 



(47) 



where H^\x) = ^^H+{x,^)^^o. The potentially diver- 
gent integral for the coefficient at the highest power of £, 
is understood as analytically regularized. Assuming that 
the forward like functions belong to the class (|52|) and 

employing the available explicit expressions for H^\x) 
one checks that the "duality property" (|47|) holds in the 
framework of the dual parametrization of GPDs. Thus all 
Mellin moments of GPD can be entirely determined from 
the knowledge of GPD in the outer region {x > £,). Now, 
e.g. employing the standard inverse Mellin transform tech- 
niques one can recover GPD in the whole region from its 
Mellin moments. In this sense the GPD in the outer region 
{x > £,) determines GPD in the whole region < a: < 1 

Let us stress once more that the possibility to fix the 
D-form factor in terms of the amplitude strongly relies on 
the postulated analyticity of the specific combination of 
the coefficients of Mellin moments of GPDs in the Mellin 
space. This analyticity assumption seems to be alluring 
from the theoretical point of view for modeling GPDs in 
the framework of the dual parametrization since it enables 
to treat spin-0 exchange contributions on the equal footing 
with that of higher spins. It also ensures the "good" an- 
alytic properties of Imyl(^) (or equivalently forward- like 
functions). However, there is still no final confidence on 
the validity of this analyticity assumption (see discussion 
in 21|). For example, this kind of analyticity can be ab- 
sent due to the contact interaction contribution to the real 
part of the DVCS amplitude or due to the so-called fixed 
pole singularity at j = — 1 (see e.g. [36]). In (|38]) these kind 
of singularities reveal themselves as terms proportional to 
a Kronecker 6-1 j which is non analytic in j. 

As it was pointed out in 14J, in the chiral quark model 
there is the explicit contact term contribution into the sin- 
glet pion GPD required by the chiral symmetry (see dia- 
gram (3) in Fig. 5 of Ref [14]) which violates the suggested 
analyticity property. 

Once the requirement of analyticity is lifted, the D- 
term may introduce an independent contribution into the 
real part of DVCS amplitude. In this case, the D-form 
factor becomes an independent physical quantity to be 
fixed from the experiment. Adding a supplementary D- 
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term 9{l — |y) SD with the Gegenbauer expansion 



SDiz) = {l~z^)Y,SdnCl{z) 



(48) 



to a GPD is equivalent to an introduction of the following 
non analytic contributions to the forward- like functions: 



x^''Q2^{x) + 26d2^-ix5{x). (49) 



Another possible source of the non-analytic contribu- 
tions to the sum rule ((55)) may be the existence of a fixed 
pole (see [551136] ) at j = —1 (i.e. angular momentum 
/ = 0). This type of contribution was reasoned with the 
Regge theory inspired argumentation for the case of for- 
ward Compton scattering amplitude in fSS' and revealed 
in the experimental measurements 39,40J. Nevertheless, 
according to 41,42,4'51Hi] via a subtracted sum rule the 
fixed pole contribution can be related to the imaginary 
part of Compton amplitude. Once generalized for the case 
of DVCS these considerations would also imply no inde- 
pendent D-term contribution into DVCS amplitude. 



4 Q2ix) and Q4(a;) from the small ^ expansion 

of Hdd{x,^) 

One of the most popular parametrizations of GPDs con- 
sists in the use of the double distributions [15,16,17,18, 
119) . In order to make this parametrization consistent with 
the polynomiality condition in its full form one has to com- 
plete it with the so-called Z?-term 0S]. More technically, 
in this case the nucleon GPD H"^ for the particular flavor 
q is parametrized in the following forrr^: 



H'^{x,^,t^0)^H%^{x,O+0{^-\x\)D'^ ^-j , (50) 

where H'^j^ is obtained as a one dimensional section of a 
two- variable double distribution F'^: 

/•I fi-\P\ 

i7?^(x,^)= dp daS{x- /3-aOF'i(P,a). 

J^i J~1+\I3\ 

(51) 

According to the proposal of Radyushkin pnirSlfTO] . the 
following model for the double distribution f is often 
used: 

F^{/3,a) = h{p,a)q{P), (52) 

where the profile function h{(3, a) is parameterized through 
the following favored Ansatz: 



r{2b + 2) [(1 



22b+lJ^2(5+ 1) (i_ 1^1)26+1 



(53) 



and q{P) {q{—l3) — — g(/3)) stands for the phenomenolog- 
ical forward quark distribution. The parameter b charac- 
terizes the strength of ^ dependence of the resulting GPD 



^ Note that since f-dependence plays no particular role in 
our considerations we just set t to zero. 



-ff'(x,^). The limiting case b ^ oo corresponds to ^ inde- 
pendent Ansatz H'^{x,£,) = 6{x)q{x) — 6{—x)q{—x), con- 
sidered e.g. in [46]. The GPD model ([STjl based on the 
Ansatz (|52|) . ([53]) for DD has been used in numerous phe- 
nomenological applications. In particular, the Radyushkin 
model is implemented in the popular VGG code (see [47]'). 

In accordance with the standard expressions for the 
singlet combination of GPDs (which is reduced to q+ {x) = 
q{x) H- q{x) in the forward limit) we define: 

H+ DD [x, = (x,0- i^x, 

pi 1-1-13 

d(3 / da {d{x - f3 - aO - S{x + (3 - a^)} 
xh('\l3,a)q+iP). 

(54) 

We now compute the coefficients at powers of ^ of A^- 
th {N- odd) Mellin moment of the quark singlet GPD 
H+dd{x,£,): 



/ dxx'^ H+dd{x,£,) 
Jo 

{N)DD , ^{N)DD^2 



= h, 

''■N+1 — " • 



,{N)DD^N-1. 
■"■N-l ? ' 



(55) 



For this we perform the straightforward calculation: 

r-l /-l 



dxx^ / d/3 / da{5{x- l3-ai)~5{x + l3-a(,)} 



1-/3 

Jo J~l+f3 

xh^'\p,a)q+{P) 

= I df3 I da^C'C'^a'=/3^-'=/j(^)(/3,«)g+(/3), 



1+13 



k=0 



(56) 



where stand for the usual binomial coefficients (C^ 



JV! 
(A'-fc)!fe! 

nl-/3 



). Now using 
daa^h^^\(3,a) 




22 r(6+iti) 



k odd; 
(1-/3)'= ,fc even; fo"" 



we obtain the following expression for the coefficients at 
powers of ^ of A^-th Mellin moment of the singlet GPD in 
the framework of the double distribution parametrization: 



(N)DD 

= I dp p^-^-'ii - prcjr^^" ' ^^^^ 



2iy+3 ^ 
2 I 



q+iP) 

(58) 



N~l 



for 1/ — 0, 1, ... 2 ■ 

Our present goal consists in deriving the expressions 
for the forward- like functions, which allow us to rewrite 
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Radyushkin double distribution Ansatz through the dual 
parametrization. For this we employ the expansion (j25|) 
of ud{x,S,) in powers of small ^ for x > Such an 
expansion was originally constructed in [17pi8j: 



H+dd{x,0 

= q+ix) 



e 



1 



92 



2! (26 + 3) dx^ 

3 



{q+{x){l^xr) 
{q+ix){l-x)^)+0{e) 



(59) 



4! (26 + 3)(26 + 5) dx^ 

<t>^oLi^) + e^fLi^) + ecp%n{x) + o(f ). 



As it was pointed out in |18j . if the singular behavior 
q+{x) '--^ ^ is assumed, the relative suppression of the 
consecutive corrections to the expansion (|59p is not 0(^^'') 
but rather 0(f ^"/x^"). As a consequence, this expansion 
is efficient only for x 3> ^. Interestingly, the expansion 
((59|) produces the correct values ([58]) of the coefficients 



and the convolution kernel Kq^{x^ y) is given by 
= -6y 



3(1-;;)- (15x + 2y) ^ 

(64) 



Note that if power-like small- a; behavior is assumed for 
q+{x): 

(65) 



the 
1 



small- a; asymptotic behavior of Q2{x) is: Q2{x) ~ 
1^. The divergency occurring in the generalized form 
factor Bio should be treated as described in sect.[3l Since 
we started from o o ([54]) with no D-term included it 
is important to check that no D-term is generated by the 
described above reparametrization procedure. 

The forward- like function Qii{x) can be expressed in a 
completely analogous way with the help of (|27|) : 



dx X 



ix) 



(N) DD 
^2iy ■ 



N -1 



(60) 



In contrast, the expression for the coefficient at the high- 
est power of ^ of the A^-th Mellin moment turns out to 
be singular. Assuming the "good" analytical properties of 
Mellin moments of H^dd in Mellin space one can treat 
this singularity exactly in the same way as described in 
sect. [3l 



dx 



„N 



(0) 



"N+lDDyX) — 



(61) 



One can check that the use of analytic regularization (|39p 
leads to the correct values /i^/^^ = 0. Thus from 
(pT|) we conclude that GPD in the framework of double 
distribution parametrization without D-term satisfies the 
"duality property" of Ref. [H]. Namely, the GPD in the 
central region < a; < ^ can be completely restored from 
its knowledge in the outer region ^ < a; < 1. 

Employing the expansions ([55]) and (|24p we establish 
the correspondence between the double distribution Ansatz 
and the dual parametrization. With the help of the general 
result (|26p the following expression for the forward-like 
function Q2ix) can be derived: 



Q2{x)^f^Q%) 



dyK^Ql{x,y)q+iy) 



(62) 



i{x) 



dyK^'l{x,y)q+iy). (66) 



The explicit expressions for fgj (x) and the convolution 

kernel KQ^^{x,y) are presented in the Appendix [Bl Un- 
fortunately, the expressions for the forward-like functions 
Q2U with v > 3 turn out to be very bulky. However, as 
it is shown in the next sect, for the case of Radyushkin 
double distribution parametrization several first forward- 
like functions provide dominant contribution into GPD 
quintessence function N{x) for small- a;. 

In order to compare the predictions of the double dis- 
tribution parametrization and the dual parametrization of 
GPDs it is extremely instructive to consider the asymp- 
totic behavior of the imaginary part of the DVCS ampli- 
tude at small values of calculated in these two models. 
Let us assume the power-like asymptotic behavior (|65p for 
q+{x). From the double distribution parametrization one 
can easily derive the leading asymptotic behavior of the 
imaginary part of the DVCS amplitude at small ^: 



ITS / ^ 



22b+i-"r(i)r(6+|)r(l- 



a) 



q+ix) 



r(2 + 26-a) 



(67) 



Let us now consider the dual parametrization. The con- 
tribution of the particular forward-like function Q21, into 
the imaginary part of the amplitude reads 



where 



/^''(-)=(2(1 



9-30a; + 29a;^ 



2 (3 + 2 6) J X 



■<l+{x) 



3- lla;^ (1 
3 + 26 



■q+{x) 



2{i-xy 

3 + 26 



j i-^i-e X 



x'^''Q2u{x) 



1 



2x 



Regge-like behavior (|65p of the forward quark distribu- 
tion leads to asymptotic behavior (|29p of the correspond- 
(63) ing forward-like function Qo for small x. For the leading 
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asymptotic behavior of its contribution into the imaginary 
part of the amplitude for ^ ~ we obtain ([28]) : 



ImA(°)(0 



2"+ir(i)r(Q + |) 

e r(a + 2) 



(69) 



Clearly, for a = b the coefficients in front of leading sin- 
gular term of ImAoDiS.) (jST)) and the Qo{x) contribution 
into the imaginary part ImA^^^ (^) coinci deQ . This fact IS 
by no doubts responsible for the approximate equality [25j 
of imaginary the parts of the singlet quark GPD contribu- 
tions into the DVCS amplitude calculated in the minimal 
dual model (with only Qo kept) to that calculated in the 
framework of the double distribution parametrization with 
6 = 1 for small values of f . Additionally taking into ac- 
count the contributions of Q2 (|62|) and Q4 (l66l) we derive 
the following asymptotic behavior for the imaginary part 
of the amplitude ImA{^): 

ImA(°)(0 +ImA(i)(e) +ImA(2)(e) 

2^+1 r{^)r{a + ^) ( (a-b) 2a{l + a) (f + a) 
^ ~^ r{a + 2) y ^ (3 + 2b) {2 + a) {3 + a) 

{a- b){a- b+ 1)2 a {1 + a) (| + a) (| + a) \ 



[3 + 2b) {5 + 2b) (4 + a)(5 + a) 



(70) 



Note, that the consideration of contributions of additional 
forward-like functions does not spoil the coincidence of 
asymptotic behavior of lmA{^) calculated from the dual 
and the double distribution parameterizations with the 
same input for b — a. Moreover, for b = a+l the contribu- 
tions of Qo and Q2 result in small-^ asymptotic behavior 
of ImA(^) coinciding with that oilmAoDiO- Analogously 
for b = a+2 one just needs to take account of the contribu- 
tions of Qo, Q2 and Q4. In general for & = a -I- M, M > 0, 
integer, it suffices to take account of a finite number of 
forward-like functions Q2U with 1/ < M obtained using 
the reparametrization procedure in order to reproduce the 
leading small-^ asymptotic behavior ([67|) of ImADD(^). 

Let us also consider the asymptotic behavior of the 
imaginary part of the amplitude for ^ ^ 1. Assuming that 



q+{x) ^ {1 — x)'* , 7>0 for 



1 



(71) 



the leading asymptotic behavior for ^ ~ 1 of the imagi- 
nary part of the amplitude calculated in the framework of 
Radyushkin double distribution parametrization is given 

by 

W..(0^2-V^ ^(^ + f;^'^ l-0^ (72) 



It is interesting to note that the factorized Ansatz with 
the correlation b — a between the profile function parameter 
b and the power characterizing the small- a; behavior of the 
forward distribution was first addressed in [18) . It was found 
to correspond to the case when the Gegenbauer moments of 
GPD are ^-independent. This very assumption was the starting 
point for the model for GPDs constructed in [48j . 



For the dual parametrization with the same input (|7ip 
the leading asymptotic behavior of ImA^"^ (^) for ^ ~ 1 is 
given by: 



ImA(°)(0 2^^ 



r(i + i 



(73) 



Contrary to the case of small-^ asymptotic behavior of 
ImA(^), for ^ ^ 1 it is not possible to make the dual and 
double distribution parametrizations with the same input 
([TT]) result in the identical asymptotic behavior of Im^(^) 
simply by adjusting the value of the parameter b. 

The two models lead to the substantially different asymp- 
totic behavior for the imaginary part of the DVCS ampli- 
tude for ^ ~ 1. In order to reproduce the asymptotic be- 
havior (|72p in the framework of the dual parametrization 
one has to take into account the contribution of the whole 
set of the forward-like functions Q2U computed with the 
help of the described above reparametrization procedure. 

Moreover, it is important to stress that the asymp- 
totic behavior ((72|) of ImA£)£)(^) for ^ ~ 1 calculated in 
the double distribution parametrization is determined by 
the value of the parameter 6, which entirely specifies the 
power-like behavior of lmA{£^) in powers of (1 — ^) for an 
arbitrary input forward quark distribution. The same pa- 
rameter b also determines the coefficient at the leading 
singular term of ImAoDiO for ^ ^ ([67|) . Thus in the 
framework of the Radyushkin double distribution Ansatz 
the two types of asymptotic behavior of ImA(^) (namely, 
for ^ ~ and <^ ~ 1 ) can not be described independently. 
This can be seen as the important drawback of this Ansatz 
making it not enough flexible for the description of the 
whole set of the DVCS data. 

On the contrary, the asymptotic behavior oilmA{£^) for 
^ ~ 1 calculated in the framework of the dual parametriza- 
tion turns out to be sensitive to the x ~ 1 behavior of the 
GPD quintessence function N{x) ioi x ^ 1: 



N{x) - 0((1 - .t)'') 



w(e)~o((i-0^)- 



(74) 



As a result, in the framework of the dual parametriza- 
tion, one can model the ^ ~ 1 asymptotic behavior of 
lmA{£^) independently from that for ~ 0. Thus, the 
dual parametrization of GPDs is more flexible and hence 
may turn to be more convenient for the description of the 
DVCS data for both large and small values of ^. 



5 Numerical results for Qo{x), Q2{x) and 
Qi{x) calculated from the double distribution 
parametrization of nucleon GPD H{x,^) 

It is extremely instructive to compare our results for the 
forward-hke functions Qo, Q2, Qi (112), (jS^l - (EH) calcu- 
lated from the double distribution model to the general 
form of GPD quintessence function N{x) (fT4|) . which can 
be recovered from the known imaginary part of the am- 
plitude Im^(^) with the help of Abel tomography method 
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1 V2x{l + x) 



2tt VTT^ 

1 x{l-x^) 



ImA 



2.T 



27r(l 



2x \3 



1 + 

2x 



ImA 



2x 



(75) 



Using the double distribution parametrization result for 
the imaginary part of the amplitude 

ImAoDiO ^jj^*' dph^"^ (^P, q+{f3) (76) 



as the input for (j75p . we easily compute the corresponding 
GPD quintessence function. 

As an example we consider the double distribution 
model for the quark singlet (C = +1) isoscalar {H\_ ') 
and isovector (i?^^'') combinations of nucleon GPDs. In 

the forward limit i/^^-* and h';^^ are reduced to the fol- 
lowing combinations of forward quark distributions: 



H^+'' {x, C = 0) = gl^^ (x) = u{x) + u{x) + d{x) + d{x) ; 

-ffj^^ (a;, ^ = 0) = q^^^ (x) = u{x) + u{x) - d{x) - d{x) . 

(77) 

As the numerical input for the isoscalar and isovector for- 
ward quark distributions q^_^''^\x) we use the LO MRST 
fit (Q2 ^ iGeV^) ;50J. 

With the help of ([62]), ^ we perform the calcu- 
lation of the isoscalar and isovector forward-like functions 
Q'^2u^^ (x) with 1/ = 0, 1, 2, which reexpress the Radyushkin 
double distribution Ansatz for GPDs in the framework of 
the dual parametrization. We compare the result to the 
general form of GPD quintessence functions N'^^'^^ {x) for 
isoscalar and isovector combinations of light quark singlet 
(C = -1-1) nucleon GPDs. It is also extremely instructive 
to compare the results for the contributions of the sev- 
eral first forward-like functions to the imaginary part of 
the amplitude lviiA^'^\S,) to the exact value of ImADD(^). 
This helps to estimate the relative importance of the non- 
forward effects encoded in Q2v{x) with v >1 for the cal- 
culation of the amplitude. 

Let us first consider the case in which the parameter b 
of the profile function (|55)) is set to its most frequent choice 
6 = 1. For this value of b one may expect that the small-^ 
asymptotic behavior of ImA£)£)(f) is sufficiently well re- 
produced already with the help of the first forward-like 



® Note that the Abel transformation was first employed in 
connection with GPDs in [TS] for constructing a model for DD 
related to PDF by Abel integral equation. More general inverse 
Radon transformation was used in [49] to relate GPDs and 
GDAs to DDs. 



function Qo since the leading small-x behavior of isoscalar 
and isovector forward-like distributions g^^'^'' ~ 1/a;"'^' ' 
is determined by the powers a'"^^ ~ 1 and < a^^^ < 1. 

In fig. [T] we show the result of an approximation of 
the isoscalar GPD quintessence function with the help of 
contributions of several first forward- like functions. We 
compare the consecutive approximations Y^=i x^'^Q2u{x) 
with W = 0,1,2 to the general result for the isoscalar 
(isovector) GPD quintessence function N'^^\x) for GPD 
in Radyushkin parametrization recovered with the help of 
Abel tomography method (|75p . One can conclude that in 
this case the contribution of several first forward-like func- 
tions really make a dominant contribution into the GPD 
quintessence functions N^^^x) for small- a;. The same con- 
clusion remains true for the case of the isovector combi- 
nation. 

In fig. [2] we compare the result for the imaginary part 
of isoscalar DVCS amplitude ImAf,^ to that calculated 

in the framework of the dual parametrization with the 

{ S) 

help of several first forward-like functions Q2J recovered 
with the help of reparametrization procedure. One may 
check that for ^ ~ the imaginary parts of isoscalar 
and isovector amplitudes calculated in the framework of 
Radyushkin double distribution parametrization of GPDs 
are with high accuracy reproduced in the framework of 
the dual parametrization of GPDs already by the contri- 
bution of the first forward-like function Qq. However, for 
f ^ 1 the contribution of several first forward-like func- 
tions Q2U into the imaginary part of the amplitude turns 
out to be insufficient to reproduce Im^^^^ with high ac- 
curacy, since the behavior of ImAf^jy for f ~ 1 essentially 

differs from that of ImAj^j^ ' ' .In order to reproduce the 
asymptotic behavior of ImA|,^ in the framework of the 
dual parametrization one has to sum up the whole series 
of contributions of Q2i'- 

Now we consider a different choice of the parameter b 
of the profile function ([53|l . We set 6 = 5 and again show 
the results for the isoscalar case. The GPD quintessence 
function A^'^'^^(a;) compared to the contribution of sev- 
eral first forward-like functions is presented in fig. [3l The 
imaginary part of the DVCS amplitude ImAf^j-, calculated 
in the framework of the double distribution parametriza- 
tion of GPDs compared to the contribution of several first 
forward-like functions into the imaginary part of the am- 
plitude is shown in fig. ID Finally, in fig.[5]we show the rela- 
tive discrepancy between ImAfyj-, and the contributions of 
several first forward-like functions into the imaginary part 
of the amplitude. In this case taking into account only Qq 
and Q2 contributions turn out to be insufficient to repro- 
duce the small-^ behavior of the imaginary part ImA^^ 
with high accuracy. In order to achieve the satisfactory 
accuracy one has to take into account the contributions of 
Qo, Q2 and Q4. 

Note, that the relative discrepancy depicted in fig. [5] 
"explodes" for f ^ 1. This is certainly due to the fact 
that the dual and double distribution parametrizations of 
GPDs result in substantially different asymptotic behav- 
ior of lmA(^) for ^ ^ 1. Hence, in order to reproduce 
accurately ImAf)]j[^) in this region in the framework of 
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the dual parametrization, one has to sum up the whole 
series of the contributions of forward-like functions that 
express the double distribution parametrization of GPDs 
through the dual parametrization. 

Finally, let us briefly consider the limiting case b = oo 
that corresponds to (2^, — 9+ i^) ■ this case in order 
to reproduce well the asymptotical behavior of ImA(^) 
for ^ ~ as well as for ^ ~ 1 it is necessary to take 
into account the whole series of contributions of forward- 
like functions Q2ij- The corresponding GPD quintessence 
function receives important contributions from Q21, with 
large v. This makes the straightforward reparametrization 
of this model for GPDs through the dual parametrization 
impracticable. 



Fig. 1. Comparison of the consecutive contributions of 
the isoscalar forward-like functions into the isoscalar GPD 
quintessence function x Qlf\x) (dashed line with short 
dashes), x {Qlf\x) + x'^Q2^\x)) (dashed line with long 
dashes), x'^ x {Q'^f\x) + x'^Q'>f \x)+x*Qi^\x)) (thin solid line) 















/ \ 




\ 








\ 





Fig. 3. Comparison of the consecutive contributions of 
the isoscalar forward-like functions into the isoscalar GPD 
quintessence function x? x Q'^^^ (x) (dashed line with short 
dashes), x^ x {Q^q\x) + x^Q'^^^ixi)) (dashed line with long 
dashes), x^ x {Q^^'^ {x)+x'^Q^i'^ (x)+x'^Q^p {x)) (thin solid line) 
into the isoscalar GPD quintessence function x? x A'^'^'(^) — 

X 111^=0^'^'' QfJi^) (tliick solid line). 































\ 








\ ^ 








\ 


/ 



Fig. 4. X ImAf,j3(^) (dotted hne) compared to the result 



into the isoscalar GPD quintessence function x^ x TV''^'*""' = for x Im^'^(^) from several first forward-like functions 



(x) (thick solid line). 



ElLoIm^'^ with W = ^ (short dashed line), W = 1 
(long dashed line) and W = 2 (solid line). 



6=; 




Fig. 2. X Imylf,j3(^) (dotted line) compared to the result 
for X Imj4'^(^) from several first forward-like functions x 
E^oIm^^'^HC) with W = (short dashed line), W = 1 
(long dashed line) and W = 2 (solid line). 



b=5 





/ / 
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/ 1 




/ / 




/ / 




/ / 













Fig. 5. Relative discrepancy between the exact result for 
lmA%]j{^) and the result for Imj4'^(^) calculated in the frame- 
work of the dual parametrization from several first forward- like 
functions ((Ml: E!1o Imyl-5('''(0 with W = 

(short dashed line), W = 1 (long dashed line) and W = 2 
(solid line). 
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6 Conclusions A On some properties of the integral 

transform Q2uix,t) H{x,^,t) 

In this paper we illustrate the application of the reparametriza- 

tion procedure, that allows to recast a particular phe- In this Appendix we consider the properties of the integral 
nomenological model for GPDs through the dual parametrizatransformation relating the singlet and nonsing let GPDf0 
tion of GPDs. We consider the popular Radyushkin dou- H±{x,£,,t) to the set of forward-like functions Q2v{y,t). 
ble distribution parametrization for GPDs and derive the 
analytical expressions for the corresponding forward-like 
functions Q2v{x), v — 1,2. It is interesting to note that 
for the most frequent choice of the parameter h of the pro- 
file function ([53|) (6 = 1) the contribution of Q2vix) with 
1/ > 3 into GPD quintessence function estimated in the nu- 
merical calculation turns out to be rather small for small- a; 
(note, that in \A\ the same property was revealed for the 
pion GPD computed in the nonlocal chiral quark model) . 
Thus, for small x the corresponding GPD quintessence 
function A''(a;) is dominated by the contributions of sev- 
eral first forward- like functions. We argue that this is due 
the fact that in this case the small-^ asymptotic behav- 
ior of ImADD(^) is well reproduced in the framework of 
the dual parametrization already with help of the few first 
forward-like functions. For larger values of h the contribu- 
tion of the forward-like functions Q2v with large v becomes 
more significant. 

We make an important observation that the forward- 
like functions Q^vix) with v >\ can contribute into the 
leading singular behavior of the imaginary part of DVCS 



The original derivation of this integral transformation was 
presented in Ref. [9j not in full details. For example, the 
precise way of treating the divergensies of the correspond- 
ing integrals has not been spelled out, which have lead to 
some confusion in the literature [5T[|28l[25ll29] . Below we 
present a more detailed and accurate derivation of this 
integral transformation. Our final result remains valid for 
the case of small-y singular behavior of the forward-like 
functions Q^viy^^^ suggested by the analysis presented in 
Sects. HIH 

First we consider the case of the singlet GPD (2;, 
that is given in the framework of the dual parametrization 
by the series: 



00 n+1 

Tl=l 1=0 



1 



1 - 



Pi 



1 
I 

(Al) 



amplitude. This provides an opportunity for a more fiexi- ^he formal series ^ corresponds to the analytic con- 
ble GPD modelling in the framework of the dual parametriza-t^';^^^^*^^^^.! 
tion. 



We also consider a way to handle divergencies, which 
may occur in generalized form-factors B^v-i o(0 due to 
the singularities of forward-like functions. The key to this 
problem is provided by the consideration of analytic prop- 
erties of Mellin moments of GPDs in Mellin space. Once 
the analyticity of the Mellin moments of GPDs is assumed, 
the so-called analytic regularization of divergencies is a 
natural way to treat the problematical form factors. Such 
approach allows to fix unambiguously the contribution of 
the D-form factor into the real part of the DVCS ampli- 
tude in terms of GPD quintessence function N{x,t) and 
the forward- like function QQ{x,t). On the other hand if 
these analyticity requirements are turned down (e.g. by 
assuming the possible fixed pole contribution or explicit 
contact term contribution) the value of the D-form factor 
turns out to be an independent physical quantity that is 
to be fixed from the experiment. 
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for the generalized distribution amplitudes entering the 
description of 7*7 hh hard process to the cross chan- 
nel (see discussion in ^23]). The crucial point is that in 
the region < ^ < 1 the representation (lAip for GPD is 
to be understood as the one which formally satisfies the 
polynomially condition {N = 1, 3, 5, ...): 



dxx'" H+{x,C,t) 



(N) 



where the coefficients h!'^\t) (with v < ^^^3-^-) at powers 
of ^ are given by the finite sums 



(A2) 



N 



2V-1-N 



(A3) 



(n 



2v+l- 



2 

{n 



+ 2) r{N + i) 



N+n ) 
2 I 



The general idea of the method 9 employed for the 
summation of the formal partial wave expansion for GPD 



For simplicity here we consider the case of spin-0 target, 
gramme of the Alexander von Humboldt Foundation, DeutschEhe corresponding formulae equally apply to the case of elec- 
Forschungsgemeinschaft. The work of K.S. was supported trie combinations of the singlet (nonsinglet) nucleon GPDs 
by STIBET scholarship from DAAD. Hf\x, t). 
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consists in presenting GPD as a result of convolution of a 
certain kernel with the set of forward-like functions Q2v{y, t) 
(y = 0, 1, ...) whose Mellin moments generate the gener- 
alized form factors Bni{t): 

Bnn+i^2u{t)^ I dyy''Q2,{y,t). (A4) 
Jo 

The explicit construction of this convolution kernel allows 
one to derive the rigorous expressions for GPDs in the 
framework of the dual parametrization. 

We start with the definition of the discontinuity of the 
particular function f{z): 

disc,=, f{z) - {fix ~ iO) - fix + ^0)) (A5) 
The basic relation derived in [5] reads 



disc 



d_ 

dy 



ds z\ 



e 1 



where 



zs = 2- 



z — 



il-s^)y' 
with < y < 1. 

Let us now consider the following function 

1 a ^ '•1 



(A6) 
(A7) 



1 



y \ dyj 7-1 ^z,2 - 2zs + e 

(A8) 

The discontinuity of this function is a suitable building 
block for the convolution kernel whose convolution with 
the set of forward-like functions Q21, allows to reproduce 
the formal series (|Aip . Let us explicitly compute the dis- 
continuity of this function employing the well known prop- 
erty of the generating function of the system of Legendre 
polynomials: 

1 1 1 



- 2Z« + Z^ 



1=0 



1 



i+i 



(A9) 



Employing (jA6[) we obtain 

disc,=,F^^-'\z,y) 



(AlO) 



The asterisk in the sum in the last line of (jAlOp denotes 
that for v = the term with ly — I = is actually absent. 
Let us consider the following integral convolution: 

V / dy {disc,=,F(2-)(z,2;) 

-disc,=^^F^^''\z,y)}Q2Ay,t)-i^^l) 



Now using (IA10|) together with the expressions for the 
generalized form factors (|A4p we obtain 



00 „i 

J2 / dy {disc.^,F^^''\z,y) 
''\z,y)]Q2uiy,t) 



j,(2v 



u=0 1=0 



00 n+1 ^ ^ , s 



(A12) 



where in the last line we have interchanged the order of 
summation introducing the new summation index n = 
2v + l-l. 

The trick is that we can compute the discontinuity of 
the function F'^'^^\z^y) with the help of an alternative 
method. Namely, instead of using the formal expansion 
(jA9|) we can consider the contribution into discontinuity 
F'^'^'^\z,y) stemming from the cut 1 — yl — < Zs < 
1 + Vl — and from the poles at Zs = for v >1. 

Let start with specifying the contribution of the cut. 
According to the standard definition of the discontinuity 
of a real analytic function (/*(z) = fiz*)) with a branch 
cut along the real axis: 



disc,=,/(z) = — ifix - lO) -f*ix- lO)) 
Zm 

^-lmfix~iO) (A13) 



Thus, the discontinuity of F'^'^^\z^ y) can be computed as 

(l\SCz2^2z,+£,^=x1-2x,+£,^F^'^''\z,y) 

\y\ oyj J-1 - 2xs + - ie \ 



J pole 
[ contrib . 



(A14) 
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Special attention is to be paid to the choice of the physi- 
cal branch of the square root in (|Al4p . According to our 
convention 

Tm ^ = Im^^ e'^ = — ^ 9{-x) (A15) 
\Jx~it ^\x\ ^\x\ 

Now let us discuss in details the contribution to the dis- 
continuity of F^^'^\z,y) resulting from poles at Zs = 0. 
Employing the expansion 



(A16) 



1 ^2xs -xi-e 



9{2x, -xi- e) 



t) 



(A18) 

As it was noted in [lOj, the sign in front of the last term 
in (jAlSp actually differs from that stated in the original 
paper [5]. 

In order to proceed further one needs to perform the 
analysis of the solutions of the algebraic equation 



2xs-xl-^^ = 0, 



(A19) 



where ; 



. The four roots of the equation (jAlOp 



Si, {i = 1, ... 4) are given by the following expressions: 



52 = ^ v/(T^^yni+M^y^TT^) 

53 = i (a - Va^^y) (l + A2)-(l-y2)) 



(A20) 



Here we have employed the notations 



1- 



The solutions of the equation si are given by y = yo 
and y — while y — yi and y = ^ are the solutions of 
the equation S3 — S4. The expressions for ?/o,i read as 



we can present F'-^'^'' {z , y) as follows: 

Obviously, for a given 1/ only a finite number of pole terms 
arise with I < 2i/ — 1. Note, that our choice of ie prescrip- 
tion is here matched with our convention (|A15P for the 
physical branch of the square root in (jA14p . Thus, we 
conclude that (|A14p may be rewritten as 

/ dyQ2Ay,t)disc,=,F'^^-''>{z,y) 

= — / dyQ2i.{y,t)- [1 + y— 
Jo V \ ay 



yo 



X (1 + fl'^ 



(1 + A^) 



'^!(i±M!)!_^2. 



^x^l + x^y 



2 V 4 

This allows to rewrite (|A18[) as 

1 

dy Q2u{y, t) (lifiCz=xF^'^''^ (z, y) 



- A2. 



(A21) 



(A22) 



^2v 
1 



{x>0 



ds- 



dy 



Vo » 
\-1v 



l-y^)g2.(2/,i) 



xi-e 



-B{~i<x<i)\ dyQ2,{y,t)-[\^y 

h y \ 



d_ 

dy 



ds- 



^2x, -x^-e 



^ ^ ^ ^ ;=o 



i-2\ ^jPiy^j 



(A23) 



A special attention is to be payed to the convergency 
of the overall integral in y in the second term of (jA23p . For 
this we need to consider the small y asymptotic behavior 
of the the elliptic integral 



^l-2iy 



ds- 



^2 



xi-e 



(A24) 



For X E {—£,'■ the following asymptotic behavior of the 
integral (jA24[l for ?/ ~ can be established: 



^l-2iy 



ds- 



y/2a 



1 

mill (-1.21^-2) 

E 

i=-2 



2iy-2 



i\ y 



'-2 \ ^jPy^j 2v-l 



X 

^2v I ^ ^ 72 



(A25) 
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where P_„(x) = ^ri-i(x)- The first term in (|A25p is ex- 
actly cancelled by the pole contribution in (jA23p . The 
term with / = — 1 in the second sum in (jA25P that is 
0{y^^) contributes solely to the D-tcrm. Finally, the term 
with I = -2 in the second sum in (|A25P that is 0(2/^"+^) 
is even in x and therefore does not survive in the singlet 
combination. 

We also have to specify our assumptions concerning 
the small y singular behavior of forward-like functions 
Q2u(,y,t)- The singular behavior of the singlet forward- 
like function Qo(y,t) is determined by that of the singlet 
combination of the corresponding forward quark distribu- 
tions: 



Qoiy) 



with 1 < a < 2 . 



(A26) 



We argue that the reasonable singular behavior of Q21, {y) 
with ly > for small y is given by 



2^(y) 



1 



with 1 < a < 2 . 



(A27) 



The next step is to add and subtract the combination 



, „ , min (-1,21^-2) , , , , 2iy-;-1 



where the functions H^\x,^^t) defined for — ^ < a; < 1 
are given by the following integral transformations: 



^ Jyo V 



l-y^]Q2.{y,t) 



rS2 , 
X / ds 

' si \/2Xs 



2 „ ^2 

1 



xt ~ 



Jo y 



i-y^ ) Q2Ay,t) 




TT X^ 



'^'^2xs-xi-e el' e 



(A31) 



with P-n{x) = Pn-i{x)- Note that the integral over y in 
the second term of (jA31l) is well convergent under the as- 
sumptions (|A26p . (jA27[) (see that the terms in braces in 
the third line of (jA3ip behave as y^"^^ for y '--^ 0) . We also 
stress that in this case the terms outside the integral stem- 
ming from (|M9l) advo cated in [5Tll25[|29| happily vanish. 
The only part of (|A30P th at sti ll may suffer from divergen- 
sies is the second sum in (jA30p containing the generalized 
form factors B2U-1 o(0- 



(A28) 



^2 



o(t) 



^ dy 



y 



y"'Q2Ay.t). 



(A32) 



from the elliptic integral in the second term of (|A23p . Now 
integrating by parts according to 



1 d 

dyQ{y)-^(ii(.y) 



dy 



dy 

d_ 

dy 



Q{y) 



m + Q{y)<t>{y) 



(A29) 



An important observation is that the second sum in (|A30[) 
is a pure Z3-term contribution. As it is explained in sect.[3l 
in order to ensure the convergency of integrals in (jA32p . 
in this case it turns out necessary to introduce a regular- 
ization (see discussion in sect. [3]). 

Our final expression (fX30l) for the singlet GPD through 
the set of forward-like functions differs from that presented 
in literature [9lll0|. The reason for this is that the results 
fOiilOj were derived under the assumption that 



and employing the asymptotic conditions (|A26p . (|A27|) we 
can rewrite the general expression (jAlip for the singlet 
GPD in terms of the forward-like functions as follows: 



limy Qo{y,t) 



0; 



lim V 



2v, 



h.{y,t)^o, iy>0. 



(A33) 



H+{x,tt) (disc,=,F(2-)(z,y) 
^=0-^-1 ^ 

- disc„F(2-)(z,y)}Q2.(y,t) 

00 

= T.^"' [Hi^\^,tt)-Hi^\-x,tt) 

+ 2 E ^ (1 - fj) (1 - fj) 4-1 (f ) S2.-1 o(i) , 

(A30) 



The first condition in (jA33p is certainly respected in our 
case, while according to the analysis presented in sects. O 
2] the second one seems to be too restrictive. However, it 
is straightforward to check that under the assumptions 
(|M3l) the resuh ^MO\i . (IMTj) for H+{x, S,, t) is reduced to 
that presented in ^lOj . 

We also present the summary of formulae for the case 
ofthenonsinglet (C = -1) GPD if_ (a;, i) ee iJ«(a;,^,t)+ 
H'^{—x,^,t). In the limit ^ ^ H-{x,£,,t) is reduced to 
q-{x, t) = q{x, t) — q{x, t). The following partial wave ex- 
pansion for H^{x, ^, t) can be written in the framework of 
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the dual parametrization: 



(A34) 



In order to sum the formal series (|A34p we introduce the 
set of the nonsinglet forward- like functions Q2v{y, t) whose 
Mellin moments give the generalized form factors Bni{t) 
analogously to (|A4[) . The reasonable singular behavior of 
the nonsinglet forward- like functions Q2v{y) with v > {) 
for small y is given by 



1 



with < a- < 1 . 



y 



(A35) 



Then, the following integral transform relating the nonsin- 
glet GPD ^, t) to the set of the nonsinglet forward- 
like functions Q2v{y) can be established: 

oo 1 

(x, e, t) = 51 / dy |disc,=,F(2'') (z, y) 
+disc,=_,F(2-)(z,y)}Q2,(y,t) 



u=0 



(A36) 



where (a;, ^, t) defined for — ^ < a; < 1 is given by 



1 - ) Q2u{y,t) 



X / ds ■ 



7^ Jo y 



1 - ) Q2uiy,t) 



„l-2iy 



withP_„(x) = P„-i(x). 



(A37) 



B Explicit expression for Q4{x) corresponding 
to the double distribution parametrization of 

GPD H+{x,0 

In this Appendix we present the explicit expression for 
the forward- like function Qi{x) calculated from matching 



of small ^ expansion of GPD H-^- in the framework of the 
dual parametrization and (j59p : 

Q^{x) = fj^l {x) + f dy K^^l {x, y) q+{y), (B38) 
where 

(6) _ / -525 + 2660 X - 5310 x^ + 5700 x^ - 2909 
^Q'^'"' " [ 2(5 + 2&) 

1365 - 4788 x + 6478 x^ - 4692 x^ + 1893 x'^ 
^ 2 (3 + 2 6) 

1 



+ (49-17.T^)(l-x^) j ^q+ix) 

-35 -h 185 X- 625 a;2 + 859x3 
5 + 26 

91 - 337 x + 625 - 603 x^ 



3 + 26 



(7-11x2) (i + x] 



5 - 70 a; + 209 a;2 17 - 94 a: + 165 a;2 



5 + 26 



3 + 26 



+ (l + a;)' 



2a;2 



-5 + 37a; 3 (3 - 11 a;) \ (1 - a;)^ (3) 



5 + 26 3 + 26 J 2x 

and the convolution kernel is given by 



(B39) 



1 



315 a; + 140 y - 270 a; ?;2 _ 84 + 27 a; y 



16 y6 
105 



(45 a; + 4 ?; — 180 xy — 16 y 
- 180 a: - 16 / + 45 a: 2/^ 
(3675 a; + 700 2; - 10500 a: ?/ 



32 (5 + 2 6) 2/*^ 
+270 a: 2/^ + 24 2/3 - 180 a: 2/^ - 16 2/^ + 45 a: 2/^ + 4 y^) 
3 



32 (3 + 2 6) y^ 
-1680 2/^ + 11010 a; 2/^ + 1344 2/^ - 5220 xy^ 
-448 2/"* + 1035 a; 2;"* + 84 2/^) . 
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